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Abstract 

A superconducting hybrid structure composed of three layers is 
considered. The 2D layers interact mutually by higher grade inter- 
layer couplings. We determine the possible superconducting modes. 
Those solutions enable to discuss the conditions for the onset and 
enhancement of 3D superconductivity in such a structure. 



1 Introduction 

We are concerned with a system of potentially superconducting 2D layers 
interacting with one another by possibly distant interlayer couplings. The 
general aim of our considerations is the study of the nanoelectronic properties 
of such structures, including the dependence of their electronic properties on 
the structural and mechanical parameters. 

We shall make use of the higher grade hybrid model (HM) of supercoductors 
formulated in [11] and developed in [12, 14, 13, 15, 16] In the present paper 
we confine ourselves to three-layer structures, with the aim to determine the 
complete set of solutions. 

According to [11], in the framework of HM the layered superconductor is 
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considered as a one-dimensional chain with "atoms" (being here identical 
atomic planes described by 2D Ginsburg-Landau theory with parameters ao 
and P) and with Josephson's bonds (called J-links) between them. The grade 
K, expressed by an arbitrary (but specified for any particular case) integer, 
defines the admitted range of Josephson's interaction in terms of interplanar 
gaps. 

We denote by ipn the order parameter associated to the layer indexed by 
the number n. Its complex conjugate (c.c.) is denoted by In the con- 
sidered variant of HM (identical atomic planes), besides the GL parameters 
characterising isolated planes there are two classes of long-range coupling 
constants: Josephson parameters 7^ and proximity efect parameters (q. In 
the present paper we intend to examine the infiuence of these long-range 
coupling, under assumption that there is no magnetic field and no currents. 

Following [11] we shall now briefly present the equations for plane-uniform 
states of our hybrid model in the absence of magnetic fleld. The order param- 
eter is then independent of the in-plane variables and the net supercurrents 
vanish. We shall confine our attention to the grade K — 2. For K — 2 the 
condition of vanishing Josephson current is equivalent to 

llii'ni'n+l - C.C.) + l2{lpn'4^n+2 + ^„-l^„+l " C.C.) = 0, (1) 

and the field equations take the form 

^ni^n + P\i'n\^i'n - ^[7l(^n+l + ^n-l) + 72(^n+2 + ^/'n-2)] = 0, (2) 

where instead of cto we have introduced 

an = ao + ^ XI C?' (3) 

with the summation runing over the planes which are J-linked to the plane 
number n. The plane index n belongs to a given set P. 

2 Exact solutions for N=3 

Let us consider the simplest non-trivial instance of higher grade system com- 
posed of a finite number of layers: the case K = 2 (hence interactions 
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of nearest and next nearest neighbours only, characterized by parameters 
Ci) 7i) C2, 72), the number of layers equal 3. We shall index the planes with 
the integers from the set P — {— 1,0,1}. The inter layer gaps will be then 
numbered by —1/2 and 1/2. The energy functional takes the form 

£ = + + aol^oP - ^[(^-1 + V'OV'o + V'o(V'-i + ^1)] 

-f (V^-iV'i) + V^iV'-i) + f (IV'il' + \W + IV'-il') 

(4) 

and the system of equations for the order parameters V^.i, and ^1 reads 

+ /3|V'ir)V'i - \bi^G + 72V'-i] = 0, (5) 
(ao + /9|V'or)V'o-^7i(V'i + V'-i) = 0, (6) 

(a_i + - ^[71^0 + 72^1] = 0. (7) 

Due to the finitenes of the system of layers under consideration, the condi- 
tions (1) for vanishing Josephson's currents J_i and Ji are automatically 
satisfied. For the present system they take the form 

7iV'_iV'o + J2ip-iipi - c.c. = 0, (8) 

7iV'oV'i + 72'0-iV'i - c.c. = 0. (9) 

Note that, according to (3), ai and q;_i are equal to one another; hence 
we shall use the symbol ai for both the cases. We have 

ai = ao + ^(Ci + C2) = CKo - ^S, (10) 

where 

5 = Ci - C2. (11) 

For the central layer we have 

ao^ao + Ci- (12) 

The trivial solution t/jj — 0, j & P describes the normal state. We are 
interested in discussing the stability of the normal solution, and in finding 
stable nontrivial solutions to the system (5-7) describing the superconducting 
states. 
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3 Properties 

The eqns. (5-7) have a number of specific properties which facihtate complete 
solution of the system. To compactify the phrasing, we shall use the following 
terminology. 

Definition 1. A solution to system (5-7) is called T-invariant iff it is 
gauge-equivalent to a real solution ipj-. 

^3 = (13) 

Definition 2. A solution to system (5-7) is called TP-invariant iff it can 
be gauge-transformed to the form 

It is also convenient to introduce the symbols for real and imaginary parts 
of the order parameters: 

V'i = % + 'i-^ji j ^ P- (15) 

Property 1. Any solution to eqns. (5-7) is gauge-equivalent to a solution 
which satisfies the following conditions 

ipQ = oo (16) 

and 

h + fc-i = 0. (17) 
Moreover, any such solution fulfills 

&i[7i«o + 72(01 + o-i)] = 0. (18) 

Property 2. Any real solution satisfies 

(ai — a_i)[Q;^ + (3{a^ -\- aia_i -|- a^i)] = 0. (19) 

where 

«t = «o + ^(72-5). (20) 
Property 3. For any solution ■j/'q = iff -01 = i^-i- 

Property 4. Any solution which is not T-invariant, is TP-invariant, and 
vice versa. 
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4 Classification of solutions 



Due to the Property 1 we can confine further consideration to the nontrivial 

solutions satisfying (16) and (17). From the remaining properties it follows 
that this set of solutions to the system (8-7) is partitoned into the following 
four disjoint classes: 



Class (A): bi = 0, ai — — a_i 7^ 0, ao — 0, 
Class (B): 61 7^ 0, ai = a_i, Oq 7^ 0, 
Class (CI): bi = 0, 0,1 = a_i, 0,07^0, 
Class (C2): 61 = 0, |ai|7^|a_i|, Oq 7^ 0. 



The classes (A) and (C2) are T-invariant, the class (B) is TP-invariant, 
and the class (CI) is T- and TP-invariant. The classes (B) and (C2) are 3D, 
the class (A) is ID and the class (CI) is 2D. In the respective classes, the 
field equations (8-7) reduce to simplified forms. Only equations of the class 
(B) contain imaginary parts. 

In the class (A) the situation is particularly simple. We have 

^ -f . (21) 

SO that (taking into account the positiveness of /3) the class is non-empty iff 

2ao + 72 - < 0. (22) 
In the remaining classes the situation is more complicated. 



Class (B): The equations take the form 



«o = -4(«o + ^^), (23) 



27: 
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where 



7i 

oi = a_i = -— Qo, (24) 
272 

bl ^ - al ^ b\, (25) 



1^1^ = 1^-1^ = -^ (26) 
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and «* is given by the eqn. (20). The necessary and sufficient condition 
for the existence of solutions in this class is the conjunction of (22) and the 
following two inequalities 

cio + < 0, (27) 

i^l - A^i)ao - 272^(72 - 6) + ^ > 0. (28) 

272 

Class (CI): The system of equations for oq and oi can be written as 

2 

ao = — (at -72 + /?a?)ai, (29) 
7i 

ai = — (ao + Pal)ao, (30) 
71 

where is given by the eqn. (20). Let us multiply by sides the equations 
(29) and (30). The result is 

{al - 72 + Pal) {do + Pal) = ^7^. (31) 
Now we make the ansatz 

do + Pal = X, (32) 

and 

2(aI-72 + /3a?)=7?-. (33) 

We are looking for non-zero solutions of the equations (29) and (30), hence 
we can divide the equations by sides. We obtain 

Pal ^ X - do> (34) 

and 

Pal=^f--al-^2>0, (35) 

X 

where both the inequalities formulate the necessary conditions for the exis- 
tence of solutions in the class (CI). Further 

iii^r - ^^ (36) 
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and, taking into account (20), 
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7i X - ao 

Geometrically, the set of solutions from class (CI) are determined by the 
points of intersection of two curves: 

y = (38) 
7i 

and 

X — Olo 

Let us note, that x — ao corresponds to the zero-solution, describing the 
transition to the normal state. We have then the equation for ao 

2(ao -S- 72)ao - 7i = 0, (40) 
which always has two real roots. We shall denote them by 

«oi ^\{S + l2- V(^ + 72)' + 87?), (41) 

and 

«02 = \{S + 72 + ^(^ + 72)^ + 87?), (42) 

For a given material, the number of solutions in the class (CI) depends on 
the temperature parameter ao in the following manner. 

1) When do < aoi, there are two solutions 

2) When ctoi < do < |((^ + 72) or |(5 + 72) < do < q;o2, 
there exists one solution 

3) For Q!o2 < CKo there are no solutions. 



Class (C2): The solutions of this class fulfill the condition 

a? + aia_i + a!.i = (43) 
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By combining it with real parts of eqns. (5)- (7) we obtain 

ai = ^(a+ + a_), a_i = ^(a+ - a_), (44) 
where the variables a+, a_ satisfy the equations 

= -4^ - 3al, (45) 

ao^--{al + ^ + (3al)a+, (46) 
2 

a+ = — (ao + /3ao)ao. (47) 
7i 

Because of (43), the inequality (22) is the necessary condition for the exis- 
tence of solutions in the class (C2). To simplify the discussion of the sufficient 
conditions let us first introduce some additional symbols: 

tti = , X = a,, K = tti (48) 

7i 

and once more change variables. We shall express the quotient and the 
product of a+ and ao by ^ and 77 defined as follows. 

> «+ 2/3 1 

C = — «i, V = ^a+oo- (49) 

ao 7i af 

Then the equations (46) and (47) will be transformed into 

e - 4«e - A77 = (50) 

and 

?7f + e + l = 0. (51) 
Eliminating rj one obtains the following equation for ^ 

e'-4^^^ + /(e + l) = 0, (52) 

which can be transformed into the form 

= {w — pf' — q, (53) 
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by introducing the variable: 



w = V6k(^ - k) (54) 

and the coefficients: 

V6X-8k^ 2 -3«;^ + «;A + A 
^=36^^'^ = ^+ 36^^ ■ ^''^ 

The request of tangency of curves 

77 = ^(e'-4<) (56) 



and 



imphes the condition 



,= -^ (57) 



2^" - - \{( + 2) = 0. (58) 

Combining (58) with (52) one can express k and A by ^: 

A = -4 — (59) 

and 

.= «5i±i. (60) 

2^ + 3 ^ ^ 

Now we can solve (60) with respect to ^: 

1 



ei,2 = -(«:-2TV(^ + 2)' + 2) (61) 

and then, apropriately substituting or ^2, calculate Ai_2 as well as r]i^2- The 

curves (56) and (57) are tangent to each other at points {^i,tii) and (^25^2)- 
For A > Ai there are no points of intersection of curves (56) and (57), hence 
no superconducting solutions to equations (46) and (47). For A = Ai there 
are two solutions. Finally, for A > Ai there are four solutions. 
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5 Onset of superconductivity 



Let us first discuss the stability of the normal state. Due to Properties the 
second variation of energy can be represented in the form 

^■'E = (ai + 2/?|^in[(<5ai)2 + {5hiY] + 

(ai + 2/3|V-ip)[(<^a_i)2 + ((56i)2]+ 

(62) 

-7i(5ai + 5a-i)5ao — 72(5ai5a_i + {5biY). 

Introducing variables a+ and a_ and their variations, we obtain for the nor- 
mal state (all fields equal zero) 

5^8 = i(ai - f )((5a+)2 + \{ai + f )((5a_)2+ 

(63) 

2(ai + l2){5biY + ao(^ao)^ - li5a+5ao. 

The necessary and sufficient condition for stability of the normal state is 
the positive definiteness of the second variation of energy, which implies the 
conjuction of the following inequalities 

ao > 0, (64) 

ao>^(5-|72|), (65) 
OiQ > q;o2, (66) 
where q;o2 is defined by the eqn. (42). 

Let us introduce the material parameters plane with coordinates 72 and 5. 
For every point of the plane (72, 8) there exists a stable normal state, de- 
pending on the temperature. The highest temperature in which (for a given 
material) the normal state becomes unstable determines the onset of super- 
conductivity. Instability of the normal state implies stability of a supercon- 
ducting state. 
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6 Stability of the superconducting states 

We shall say that a solution from a given class has the property of internal 
stability if it is stable with respect to variations preserving the class. 

For the classes with real solutions we have 

5^8 = {ai + 3(3al){6aif + (ai + 3(3al^){6a_if+ 

{ao + 3/3a2)(5ao)2 + [ai + /5(a? + al,) - 72](5(&i)'+ (67) 

— 7i(5ai + Sa^i)Sao — 725ai5a_i. 
hence for the class (A), taking into account the solution (21) 

S'^S = -{di + 72)(5a+)^ - (tti + —){da_Y + ao{6aof - ^i6a+Sao. (68) 

Note that the coefficient by ((^&i)^ equals zero. The conditions for stability 
of mode (A) read 

ao < 2 ~ 72, (69) 

ao<^(5 + 72), (70) 
b 1 

«o(<5o-(--72)) + ^7?<0, (71) 

and the condition (64). To fulfill (69) and (71) the following restriction of 
the material parameters is necessary 

2-72>|7i|- (72) 
For such materials, the real roots of the polynomial from (71) have the form 



«i,2 = ^(^-72±y(^-72)^-7n. (73) 

Within the class (CI) the second variation of energy has the form 
= (2ai + 72 + 3/?a^) {baSf + (2ai - 72 + 6/3a^) ((5a+)2+ 

(74) 

(2ai - 72 + 2/?a^)(56i)2 + {olq + 3/?a^)(5ao)2 - 27i5a+5ao. 
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and the stability conditions read 

2do - 5 + 72 + S(3al > 0, (75) 

2do - (5 - 72 + 2f3al > 0, (76) 

do + 3pal > 0, (77) 

(2do - 5 - 72 + 6pal) {do + Spal) - -ff > 0, (78) 

The conditions of positive definiteness of the second variations for the 
class (C2) 

S^S = {2di + 72 + 3(3{al + 2al)){Sa_)^ + (2ai - 72 + 3(3{2al + al)){Sa+y+ 
(2di - 72 + 2(3{al + a'i)){5hy + {do + 3(3al){5ao)^- 
27i5a+5ao + 6/3a+a_5a+5a_, 

(79) 

and for the class (B) 

S^E = (2ai + 72 + (3{3al + 2bl)){5a.f + {2d^ - 72 + 2p{3al + bl)){6a+f+ 
(2<5i - 72 + 2p{al + 3bl)){6biy + {do + 3Pal){6aof- 
2jiSa+Sao + 2j3a+biSa+dbi, 

(80) 

are necesarry and sufficient to ensure the stability in the corresponding 
classes. Those modes, however, do not appear at the onset of the super- 
conductivity in the system under consideration. 

7 Conclusion 

1. The isolated atomic planes are all described by the model 2D GL with the 
same parameters cto and /?. The parameter ao depends on the temperature. 
If ao > 0, then the 2D state of such a plane is N (normal), in the oposite 
case it is the state S (superconducting). We interpret the parameter ckq as a 
mesure of empiric temperature, and introduce the notation t — To + ckq- The 
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result for an isolated plane: state N for r > tq, state S for r < tq. 

2. Each material is charcterised by two empiric temperatures ta and tc 
defined by the formulae 

TA = To- + ^(5-72), (81) 

and 

Tc = To - Cl + «02, (82) 

where 0:02 is given by the eqn. (42). As long as r > maa;(r4, tc), the normal 
state is stable. Below this limit a stable superconducting mode apcarcs. If 
Ta > Tc, the mode is A; in the opositc case - it is the mode C. The other 
superconducting modes do not appear in the onset. 

3. On the material plane (72, 6) the limiting curve between the regions 
of onset A and onset C is placed at the half-plane 72 < 0, and given by the 
formula ^ 

5 = 72 - -• (83) 
72 

4. The solutions belonging to the class (A) exist in the whole plane (72, S). 
The unique condition for their existence is (5o < |(<^ — 72), hence r < ta- 
This condition ensures also the internal stability of the solutions. 

5. Under the above circumstances the necessary and sufficient condition 
for the stability of solutions from the class (A) is their C-stability expressed 
by the eqn. (71). However C-stability of the solutions from the class (A) is 
possible only for the materials fulfilling the inequality (72). 

6. Then the necessary condition of such stability is ao > ai, (with oi 
given by (73)), hence r > tca- Violation of the condition implies the loss of 
stability on behalf of the mode (C2). 

7. If material is of the type (A), the solutions from the class (A) are 
stable in the interval ai < ao < 0:02, hence in the interval tqa < t > ta- 

8. If the material is of the type (C), then the solutions from the class (A) 
are stable in the interval ai < cio < 0:2, hence tca < t < tac- At the both 
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ends the stability is lost on behalf of the mode (C2). 

9. Now we shall discuss the enhancement of superconductivity in the sys- 
tem under consideration. We consider the enhancement from the following 
two points of view: 

a) The enhancement due to long distance coupling with respect to the 
short distance ones, 

b) The enhancement of 3D superconductivity with respect to the 2D. 

10. Consider first increments of the onset temperatures ta and tq due to 
long distance couplings. Let taq and tcq denote the onset temperatures ta 
and Tc for the first grade material, i.e. for C2 = and 72 = 0. Then, from 
the equations (81)-(82) it follows that 

TA-TA0 = -^(C2 + 72), (84) 

and 

TC - Tco = -^[V(Ci- ((2-72^ + 871^ - VCf + 87? + (C2 - 72)]. (85) 

11. In consequence of the above formulae, the mode A superconductivity is 
enhanced provided that 

C2 + 72 < 0. (86) 

Similarly, due to the fact, that the right hand side of the eqn. (85) is a 
monotically increasing function of the difference C2 ~ 72, the mode C su- 
peconductivity is enhanced provided that 

C2 - 72 > 0. (87) 

In consequence, negative values of 72 are in favour of both the modes A and 
C. At the same time, negative values of ^2 enhance the mode A and tend to 
supress the mode C of superconductivity. 

12. For some materials, charactised by appropriate values of the parame- 
ters Ci, C2 and 72, the formulae (81) and (82) can give the onset temperatures 

Ta and/or tq greater than the 2D critical temperature tq. In such a situation 
the out-of-plane superconductivity appears in spite of the fact that all the 
layers remain in overcritical in-plane states. 
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